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Abstract 

We consider a cellular satellite system conceived to en- 
hance the capabilities of the pan-European terrestrial sys- 
tem (GSM). This adopts EHF band and highly-inclined 
orbits . We present a preliminary assessment of system 
capacity based on asynchronous direct-sequence spread- 
spectrum multiple access (DS-SSMA). Performance is mea- 
sured in terms of error probability achieved by A users 
simultaneously accessing the system with a given signal- 
to- noise ratio . 

INTRODUCTION 

A cellular satellite system has recently been pro- 
posed [3, 4, 12] which is conceived to enhance the capa- 
bilities of the pan-European terrestrial cellular system 
(GSM). Intended as an efficient complement to GSM, 
this new system adopts EHF band (specifically, 40/ 
50 GHz), and highly-inclined orbits. Its main goals 
are: [3, 12] (i) To provide higher source data rates (64 
kbit/s, and possibly n x 64 kbit/s) in order to support 
ISDN-compatible teleservices, (ii) To expedite the dif- 
fusion of mobile services in Europe and expand their 
coverage to Eastern Europe, part of North Africa, and 
the Near East, (iii) To integrate maritime and aero- 
nautical users into the terrestrial network, (iv) To offer 
private users a solution for implementation of mobile- 
radio closed networks with coverage radii up to a few 
hundred kilometers, (v) To allow implementation of 
networks of small, portable terminals in the perspec- 
tive of a forthcoming worldwide personal communica- 
tion system. 

In this paper we present a preliminary assess- 
ment of system capacity based on asynchronous 
direct-sequence spread-spectrum multiple access (DS- 
SSMA). System performance is measured in terms of 
error probability achieved by K users simultaneously 
accessing the system with a given signal-to-noise ratio. 


tern [13],' K users transmit digital data over a common 
satellite channel with center frequency f c . The signal 
transmitted by the Jc th user is 

s fc (<) = V2P b k (t)a k (t) cos(2ir f c t + 9 k ) ( 1 ) 

where P denotes the common transmitted power, 6*(<) 
is the data-bearing signal, afc(t) is a periodic waveform 
with period T and formed by N rectangular chips, each 
with duration T c ~T/N and amplitude ±1, and Ok is 
the carrier phase. The signal bk(t) is a rectangular 
waveform taking on the value bkj in the time interval 
[IT, (l + 1)T], where bkj the /th source bit from user 
k, taking on values ±1 with equal probabilities and 
independently of the other 

The code sequence assigned to the k th user takes on 
value a"p in the interval [ jT c , (j-f 1)T C ] for a code word 
a^\ a^ -1 ^ taking values ±1. The signal 

bandwidth is proportional to = N/T, and hence 
N is the bandwidth-expansion factor. 

The signal received by the kth user is 

K 

'J’lP y l a k (t - T k )b k (t - r k ) cos(27r/ c f + 4>i fc) + n k (t) 
k = 1 

where n*(J) is a white Gaussian noise process indepen- 
dent of the other random variables (RV) and whose 
two-sided power spectral density is Nq/2. The de- 
lay rjt and the phase <j>k = Ok — 2 nf c T k model asyn- 
chronous transmission. We assume that 7> is uni- 
formly distributed in [0, T], <j>k is uniformly distributed 
in [0, 27r], and r*, 4>k are independent of each other and 
of the other RVs in this model. 

The ibth-user receiver computes the correlation be- 
tween the received signal and ajt(f) cos(2^/ c f) in the 
time interval [0, T], and makes a decision on the value 
of bk, o based on the sign of this correlation. Consider 
receiver #1, and with no loss of generality assume 
T\ = <j)\ — 0. With f c 1/T, the correlator output is 


SYSTEM MODEL Zl = J^Tb lfi + J^u + n, (2) 

In the model of an asynchronous DS-SSMA sys- ^ ^ ^ ^ in tfce nght . hand side is the use _ 

ful signal, the second term is the interference from 
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the other users, and n is white Gaussian noise. The 
multiple- access interference RV u is defined as 


u = + h,oRk t i(rk)\ cos^jt = 

k-2 

K 

= X wk cos < * >k 

k-2 

(3) 

where and R k ,i are partial crosscorrelations 

Rk,i(r) = f a k (t - r)ai(t)dt 

h ( 4 ) 

Rk,i(r) = J a,(t - r)ai(t)dt 

with 0 < r < T, 

The average error probability for user #1 can be 
evaluated by using the moments of this RV u. 

Assume 6 lj0 = 1. Then the average error probability 
P(E) is equal to P(Z\ < 0), and we have 

p(£,= /°[v^( 1 + f)] / “ ( “ )<l (5) 


where PT/N 0 — Eb/No is the signal- to-Gaussian noise 

power ratio, and Eb = PT is the energy per bit. The 
function Q(x) is defined in the usual way: 


1 


where f x (x) is the pdf of the RV x. This can be ex- 
panded in the form 

| / ^,2 \ 00 

/*(*) = ~^= ex P ( - t) Yl CnHen W ( 8 ) 


where He n (x) is the nth Hermite polynomial 
He n (x) = (-!)" exp T^exp (~T 


The coefficients c n are obtained from the orthogonality 
relation of Hermite polynomials, and are given by c n = 
E n /n\ , where E n = E[He n (x)\. The coefficients E n 
are generalized moments, and can be written in terms 
of the central moments of x. Since the pdf of x is an 
even function, and He n (x ) is odd when n is odd, E n 
is zero for n odd. The generalized moments Z? 2 n can 
be expressed in the form 

n 

E 2 u = X h ^, 2 kE[x 2k ] 

n = 0 

where /i2n t 2fcj the coefficients of the orthogonal poly- 
nomials Hc 2 n (*^)i cari he computed using the recursion 

L n(2n - 1) L , n , 

^2n,2Jfc = 7 ^2n — 2,2fc , » = U, . . . , 71 — 1 


/*2n,2n — T 


for all n 


The coefficients jE^n can be computed from [13] 


E[x 2k } = zjj 


( 2JA /£[u 2 *] (2 k — 2 t)! / Nq 
^ \2i) T 2i 2 k ' i (k-i)\ \2PT 


Once the moments of the RV Z\ have been eval- 
uated, several techniques are available to compute 
the error probability (5). Refs. [1, 2, 14] summarize 
some of them, based on series expansions or on Gauss 
quadrature rules. In the following we shall use a tech- 
nique advocated in [13] and based on Gram-Charlier 


By substituting (8) into (7), and using the previous 
results, we obtain a series expansion for the average 
error probability: 

P(E) = Q(*o) - e-*o / 2 X TTTTTT 

n= i (2nj! v27t 


Gram-Charlier expansion 

Gram-Charlier series expand a probability density 
function (pdf) in terms of derivatives of a known 
pdf [5], 

Following [13], define a RV x associated with 
multiple-access interference and Gaussian noise: 


E[u 2 ) No 
T 2 + 2 PT 


x = xo — h 
V u 


where Xq is a normalization factor for x. The average 
error probability is given by 


P{E) = P{x < -* 0 ) 


(x)dx (7) 


The first term is the error probability we would 
achieve in the presence of an additive Gaussian dis- 
turbance whose power spectral density is the sum of 
the multiple-access interference and of the Gaussian 
noise. 

This technique provides satisfactory results when 
the number of users is large and the signal-to-Gaussian 
noise ratio is low [13]. In this situation, two to three 
terms in the summation are enough for a good approx- 
imation to error probability. The moments of the RV 
u can be computed as described in [14]. 

Gaussian assumption 

When the number of users is large enough, the pdf of 
the multiple interference is often assumed to approx- 
imately Gaussian, with a variance equal to the sum 


462 



of the variances of the multiple-access interferes and 
of the Gaussian noise. The average error probability 
based on the Gaussian assumption can be expressed 
as 


P{E) = Q 


( No K- I V* 
\2E b 3 N ) 


( 11 ) 


Choosing the code sequences 

A class of binary code sequences with good correla- 
tion properties is provided by the maximum-length se- 
quences. Another choice is suggested by Pursley and 
Eoefs [10]: among all the maximum-length sequences 
with period N = 2" - 1, we retain only those whose 
correlation sidelobes have lower energies: these se- 

quences are called AO/LSE (auto-optimal/least side- 
lobe energy). The main drawback with these se- 
quences is the small number of them available for each 
value of N . Another possible choice is provided by 
Gold sequences [7]. 


NUMERICAL RESULTS 

In this work we consider a pseudo-random sequence 
with length N = 255 generated by an 8-stage shift 
register whose feedback connections are described by 
the polynomial 453 (octal notation). We assume that 
each user is assigned a shifted version of the sequence, 
so that a maximum of 255 users can be accommodated 
by the system. 

Fig. 1 shows the average error probability of the 
system versus the number of users for several values 
of signal-to-Gaussian noise ratio E b /No. The curves 
are obtained through the Gram-Charlier expansion 
method described before. The asterisks denote the ap- 
proximation obtained under the Gaussian assumption, 
and refer to the error probability curve lying above 
them. It can be seen that the Gaussian assumption 
provides accurate enough results when the number of 
users is very large, and hence for high P(E ) values. 
The curves were obtained by retaining three terms in 
summation (10), and taking moments of the RV u up 
to order six. 

Fig. 2 reorganizes the results by showing the number 
of users that can access the system at the same time, 
as a function of the signal-to-Gaussian noise ratio for 
several values of error probability. We can see that, if 
an error probability less that 10” 6 is sought, we need 
an Eh l No of at least 20 dB for 25 simultaneous users. 
If each user transmits at a speed of 64000 bit/s and 
the processing gain is 255 (24 dB), the bandwidth used 
is 64000 x 255 = 16.3 MHz, that is, each user needs 
a bandwidth of about 650 KHz. Under these condi- 
tions CDMA multiple access offers no advantage with 



Figure 1: Average error probability P(E) versus the 
number of users for different values of the ratio E b /No 
for maximum-length sequences with N = 255. 



Figure 2: Number of simultaneous system users vs. 
E b /N 0 . The length of the code sequences is 255. 


respect to FDMA, which would need about 128 KHz 
to transmit at 64 Kbit/s. Thus, SSMA would need 
an efficient coding scheme (which provides a suitable 
coding gain) to become attractive. 

To increase the number of users that can access the 
system while keeping the same error performance we 
may think of using longer sequences, which will in- 
crease processing gain. The average error probability 
P(E) obtained with a sequence of length 511 (feedback 
shift register with nine stages and generating polyno- 


463 





mial 1021) was computed. The corresponding values 
of system capacity vs. E b /N Q are shown in Fig. 3. We 
see that to achieve K = 25 and P(E) = 10~ 6 as be- 
fore, 13 dB are enough, with a savings of about 7 dB. 
The cost is a twofold increase in bandwidth occupancy. 



Figure 3: Number of users that can access the sys- 
tem with an assigned performance vs. E b /N o- The 
sequence length is 511. 


Consideration of the satellite channel 

We now consider some of the major factors that in- 
crease the capacity of an SSMA system [6], viz., 
voice activity, spatial discrimination provided by cross- 
polarization frequency reuse, and multi-beam anten- 
nas. 

Voice activity. With a 35% voice activity factor, 
the number of users in Fig. 2 and Fig. 3 must be mul- 
tiplied by 1/0.35 = 2.86. Thus, when N = 255, for 
P(E) — 10 -6 and E b /No = 20 dB we have about 70 
users, and the bandwidth is 650000 x 0.35 = 227.5 Kz 
per channel. With N = 511, for the same values of 
P{E) and E b /No we have about 140 users, and the 
bandwidth per channel is about the same. 

Cross-polarization frequency reuse. As for fre- 
quency reuse by cross-polarization, a quantitative 
analysis can be performed by using the computational 
techniques described before. Assume that K users ac- 
cess the system with each polarization. If X p denotes 
the relative signal level with a given polarization, and 
received by a user transmitting with the opposite po- 
larization, the total signal received by the fcth user is 

K 2 K 

r l (() = ^- 4 *( < ) + A P^’ E Mt). 

k = \ k-K+1 


where 

Ak(t) = dk(i - T k )b k {t - Tjt) cos(2:r f c t + <j> k ). 

The output of the fcth-user correlator is 

Z i = ]/j Tb i,o+ \fj u + Xp\f^u' + n (12) 

where the first term in the right-hand side is useful 
signal, the second the interference from the users with 
the same polarization, and the third the interference 
from the cross-polarized users. As usual, n is Gaussian 
noise. The RV u was defined in (3), while u ' is defined 
as follows: 

2 K 

«■ = e h + b k,oRk,l(Tk)] COS 4> k 

k-K + 1 

2K 

= E U ’ k COS 

k = K + l 

(13) 

The two RVs u and u 1 are independent, because b] C) 
rjk, and <f>k are independent for each user. If we de- 
fine the new RV v — u -{- X p u\ we obtain for error 
probability the same expression as (5), and hence to 
compute it we can use the same technique as before. 
To evaluate the moments of t; we follow [8], and obtain 

E[v 2n ] = it (ll) E[u 2p ]X 2n ~ 2p E\u' 2n ~ 2p '\ (14) 

P — 0 ' 

where the moments of u and u* are calculated as above. 

The error probability obtained with a cross- 
polarized channel attenuation equal to 6 dB, and se- 
quence length 255 was computed. The corresponding 
values of capacity vs. E b /No are shown in Fig. 4. We 
can see that the capacity increase is about 85%, which 
could not be obtained by FDMA. 

Fig. 5 shows capacity curves with sequence length 
255, cross-polarization diversity, and a 35% voice ac- 
tivity factor. 

If error-control coding is used, these curves move 
toward left by an amount equal to the coding gain. 

Spatial discrimination Channel capacity can be 
further increased if multi-beam antennas are used, 
which allow a degree of frequency reuse. The amount 
of reuse depends on the maximum interference that 
can be tolerated among beams transmitting in the 
same bandwidth. A factor of frequency reuse equal 
to 2 or 3 seems reasonable, which implies increasing 
the capacity by a factor 2 or 3 if the interference is 
negligible. 
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Figure 4: Number of users that can access the system 
with a given P(E) for an attenuation of the cross- 
polarized channel equal to 6 dB. 



Figure 6: System capacity per beam for relevant values 
of P(E) and E b /N 0 . 

number of channels. By doubling the bandwidth, the 
achievable improvement is only 0.5 dB. 



Figure 5: System capacity with single beam for several 
values of average error probability, an attenuation of 
the cross-polarized channel equal to 6 dB, and 35% 
voice activity factor vs. E^/Nq. 

SYSTEM CAPACITY 

Consider now the satellite mobile radio system de- 
scribed in the Introduction. Assume 1000 channels 
with 40 beams, i.e., 25 channels per beam. From 
Fig. 6, obtained from Fig. 5 by zooming in the relevant 
region of E b /No and P(E) values, we can see that for 
25 users to access the system with P{E) — 10“ 6 we 
need E^/Nq ~ 11,3 dB. With code sequences of length 
N = 511, we need about 10.8 dB to achieve the same 


A note on synchronous SSMA 

To achieve a more efficient use of the bandwidth with 
spread-spectrum modulation, multiple-access noise 
due to asynchronous transmission should be elimi- 
nated. This may be accomplished by making the 
users’ spreading modulation synchronous at the satel - 
life , thus providing channels that are ideally orthog- 
onal, In synchronous design, the users’ signal would 
have to be synchronized to a small fraction of a chip 
time to achieve the desired orthogonality and mini- 
mization of interference. In this system, timing con- 
trol algorithms, their accuracy, and their impact on 
system design are be important issues that will not be 
dealt with here. In the following, capacity evaluations 
are extended to synchronous SSMA. 

With a fully synchronous system the delays r % in 
the system model can be disregarded, and we assume 
rjfc = 0 for k = 1, 2, . . . , K . The signal received by the 
k th user is 

K 

n(t) = \PlP y ] a k (t)b k (t) cos(2n f c t + 9 k ) + n k (t) 

k = 1 

(15) 

and the output from the correlator is 

Zi = \f^ b i.oT+ \Jjw + n (16) 

where the RV w is defined as 
K 

w = y>-,ofl M (0) cos 9 k (17) 

fc = 2 
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with 


Rk, i(0) = / afc(O a i(O di * ( 18 ) 

J o 

Apart from a factor 1/T, is the crosscorrelation 
between a pseudo-random sequence and its translated 
version, . It is well known that Rk t i(0) takes values 
—T/N = — T c , with N the sequence length. Thus, the 
EV w can be given the form 

K 

iv = T ^6^0 cos 0* (19) 

k-2 

and is the sum of K — 1 independent RVs. 

To evaluate error probability, we compute the mo- 
ments of w and apply the Gram-Charlier expansion 
method as before. Fig. 7 compares the error probabil- 
ity of a system with 250 users and that of a chan- 
nel with a single user and additive white Gaussian 
noise as the only disturbance. If an error probability 
10 -6 is desired, we need only to increase the signal-to- 
noise ratio by 0.2 dB with respect to the Gaussian 
channel. The spectral efficiency increases consider- 
ably, as we can allocate 250 users in the bandwidth 
64000 x 255 = 16,3 MHz. 



Eb/No [dB] 


Figure 7: Average error probability vs. E^/Nq with 
synchronous transmission and sequence length 255. 
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